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Abstract. A remarkably quantitative understanding of the electrical and mechan- 
ical properties of metal wires with a thickness on the scale of a nanometer has been 
obtained within the free-electron model using semiclassical techniques. Convergent 
trace formulas for the density of states and cohesive force of a narrow constriction 
in an electron gas, whose classical motion is either chaotic or integrable, are de- 
rived. Mode quantization in a metallic point contact or nanowire leads to universal 
oscillations in its cohesive force, whose amplitude depends only on a dimensionless 
quantum parameter describing the crossover from chaotic to integrable motion, and 
is of order InN, in agreement with experiments on gold nanowires. A linear stability 
analysis shows that the classical instability of a long wire under surface tension can 
be completely suppressed by quantum effects, leading to stable cylindrical config- 
urations whose electrical conductance is a magic number 1, 3, 5, 6,... times 2e^ /h, 
in accord with recent results on alkali metal nanowires. 



1 Introduction 

In 1971, Gutzwiller's trace formula |^ expressing the quantum density of 
states of a classically chaotic system as a Feynman sum over classical periodic 
orbits gave birth to the field of quantum chaos. In the subsequent decades, 
the trace formula was generalized, and applied to a wide variety of physical 
systems |2]j^ . Of particular interest here are trace formulas for systems with 
continuous symmetries and broken symmetries One of the most 

important successes of this semiclassical approach has been the description 
of shell effects in finite fermion systems ^,||. In this article, we discuss the 
application [p| pO| , pl] , p^ of trace formulas to describe quantum-size effects 
in a particular class of open quantum systems: metallic nanocontacts and 
nanowires. 

In the past eight years, experimental research on atomically-thin metal 



wires has burgeoned |T30,|l|,[§0,|l|,|l|,|gj2|,|2| . In a seminal experiment 
[ p!4| carried out in 1995, Rubio, Agrait and Vieira simultaneously measured 
the electrical conductance and cohesive force of an atomic-scale gold con- 
tact as it formed and ruptured (see Fig. [^). They observed steps of order 
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Fig. 1. Left: Simultaneous measurements of (a) tiie conductance and (b) tlie colie- 
sive force of a gold nanowire during elongation at room temperature, from Ref. |^| . 
Right: Transmission electron micrographs of an atomic-scale gold contact breaking, 
from Ref Q. The measured electrical conductance of the contact is (d) G ~ 2Go, 
(e) G ~ Go 



Go = 2e^//i in the conductance, which were synchronized with a sawtooth 
structure with an amplitude of order InN in the force. Similar results were 
obtained independently by Stalder and Diirig jisj. For comparison, electron 
micrographs by Ohnishi et al. pof illustrating the atomic-scale structure of 
a gold nanocontact breaking are also shown in Fig. ^ 

Conductance steps of size Go were first observed in quantum point con- 
tacts (QPCs) fabricated in semiconductor heterostructures |2^, and are a 
rather universal phenomenon in metal nanowires jlSj, even being found in 
contacts formed in liquid metals p7[ | . The precision of conductance quantiza- 
tion in metal nanocontacts is poorer than that in semiconductor QPCs due 
to their inherently rough structure on the scale of the Fermi wavelength \p, 
which causes backscattering [Q, and due to the imperfect hybridization of 
the atomic orbitals in the contact, especially for multivalent atoms [Q. As 
we shall see in the following, the sawtooth structure in the cohesive force can 
be considered a mechanical analogue of conductance quantization p5| . 

A remarkable feature of metal nanowires is the fact that they are stable 
at all. Fig. || shows electron micrographs by Kondo and Takayanagi illus- 
trating the formation of a gold nanowire. Under electron beam irradiation. 
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Fig. 2. Transmission electron 
micrographs showing the for- 
mation of a gold nanowire, 
from Ref. [Q; (a) an image 
of Au(OOl) film with closely 
spaced nanoholes, the initial 
stage of the nanowire; (b) a 
nanowire four atoms in diam- 
eter, resulting from further 
electron-beam irradiation 



the wire becomes ever thinner, until it is but four atoms in diameter. Almost 
all of the atoms are at the surface, with small coordination numbers. The 
surface energy of such a structure is enormous, yet it is observed to form 
spontaneously, and to persist almost indefinitely. Even wires one atom thick, 
such as that shown in Fig. |l|(e), are found to be stable for days at a time 
[ p^[|20| . Naively, such structures might be expected to break apart due to 
surface tension but we shall show that quantum-size effects can 

stabilize arbitrarily long nanowires fl^ . 

2 Free electron model 

We investigate the simplest possible model |2^,^ for a metal nanowire: a 
free (conduction) electron gas confined within the wire by Dirichlet boundary 
conditions. A nanowire is an open quantum system, and so is treated most 
naturally in terms of the electronic scattering matrix S. The Landauer for- 
mula ||30| , ^ expressing the electrical conductance in terms of the submatrix 
5*12 describing transmission through the wire is 

G=^/ dE^^Tr{sl,{E)SME)}, (1) 

where f{E) is the Fcrmi-Dirac distribution function. The conductance of a 
metal nanocontact was calculated exactly in this model by Torres et al. . 
The appropriate thermodynamic potential to describe the energetics of such 
an open system is the grand canonical potential J7, whose derivative with 
respect to elongation gives the cohesive force F: 

f2 = -^JdEgiE)ln{l + e-^^^-^^), F = (2) 

Here f3 is the inverse temperature, /i is the chemical potential of electrons 
injected into the nanowire from the macroscopic electrodes, and g{E) is the 
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contact diameter [Xp] 



Fig. 3. Electrical conduc- 
tance and cohesive force of a 
nanowire, modeled as a nar- 
row neck in a free-electron 
calculated from Eqs. 



(E|)~0 ^^'''^ temperature. 
The S-matrix was calcu- 
lated using the adiabatic and 
WKB approximations, fol- 
lowing Ref. [^. For compar- 
ison, the contribution to the 
force from the surface ten- 
sion and curvature energy is 
shown as a dashed line. Note 
that £f/Af = 1.7nN in gold 



electronic density of states (DOS) of the nanowire. The DOS of an open 
system may be expressed in terms of the scattering matrix as |Q 

m-^Tr[sHE)^-n.c.]. (3) 

This formula is also known as the Wigner delay. Note that in Eqs. (|l|) and (||), 
a factor of 2 for spin degeneracy has been included. Thus, once the electronic 
scattering problem for the nanowire is solved, both the conductance and force 
can be readily calculated | p5| , p9| , p4[ , as shown in Fig. |[ One sees that there 
is an almost quantitative agreement with the experimental results shown in 
Fig. |l|: for example, the force necessary to break the last conducting channel 
is approximately sf/^f (=1.7nN in gold), where ef is the Fermi energy. 

3 Weyl expansion 

In order to separate out the mesoscopic sawtooth structure in the force, 
associated with the opening of individual conductance channels, from the 
overall (macroscopic) trend of the contact to become stronger as its diameter 
increases, it is useful to perform a systematic semiclassical expansion l^jsj of 
the DOS, g{E) = g{E) + 6g{E), where 5 is a smooth average term, referred 
to as the Weyl contribution, and 5g{E) is an oscillatory term, whose average 
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is zero. For the free electron model with Dirichlct boundary conditions, the 
Weyl term is [|| 



where ks = ^/2mE/h, V is the volume of the wire, A its surface area, and K 
the integrated mean curvature of its surface. Inserting Eq. (^) into Eq. (||), 
one finds the following semiclassical expansion at zero temperature: 

n _ 2klV kj,A 2kFK 5Q 
7^^~ 157r2 " 97r2 + 77' 

One can show jl^ that interaction effects are higher order in h. In the same 
spirit, a semiclassical expansion for the conductance G = {2e^/h)Gs + SG 
may be derived, where the corrected Sharvin conductance is |32| 

Here D* is the narrowest diameter of the nanowire. 

When the wire is elongated, the atoms rearrange themselves, but the vol- 
ume per atom remains essentially constant ||l0| , p^ . Thus, when differentiating 
Eq. (H) to calculate F, the first term on the r.h.s. gives zero: 

The cohesive force is given by surface tension plus a small curvature correction 
(the sum of which is indicated by a dashed curve in Fig. ^), combined with 
an oscillatory quantum term. 



4 Trace formulas 

The oscillatory contribution Sg{E) to the DOS may be approximated as a 
Feynman sum over classical periodic orbits a la Gutzwiller ||^,^ . Since we are 
interested in modeling nanowires which may possess axial and/or transla- 
tional symmetries, however, we can not in general utilize Gutzwiller's original 
trace formula jl], which describes systems whose periodic orbits are isolated, 
but must instead employ a generalization due to Creagh and Littlejohn, de- 
scribing a system with an /-dimensional Abelian symmetry 1^: 

r /^x 2 1 ^ TrVrJr^'^ f Sr 

^■'^(^^ - ^(2^^ IdetM. - 11 V. U - — - Tj ' 



where the sum runs over /-dimensional families F of degenerate periodic or- 
bits, Tr is the period of an orbit in T, Vr is the /-dimensional volume spanned 
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by r, Sr is the action of the orbit, and cr/^ is a phase shift determined by the 
singular points along the classical trajectory. The quantity M is the so-called 
monodromy matrix, characterizing the stability of the orbit with respect to 
perturbations. It describes as a Poincare map the linearized motion of small 
perturbations from the periodic orbit in a surface of section perpendicular to 
the orbit in phase space: an initial variation of momentum and position in 
the surface of section {Sr,Sp) is related to the mismatch {6r',dp') after one 
period by 



Sp' J \6p 

Finally, the factor Jp — \ det{dr' /dp)\. 

We shall also need to consider the breaking of continuous symmetries, 
which is elegantly described in terms of semiclassical perturbation theory 
wherein the cosine in the trace formula is replaced by 



cosiSr/h + Or) ^ Re {e'(^^/'*+^^) (e'^^^^")^} ' (10) 
where 

^^^^ASr/r^^^ ^ y-1 J d^(^gyASr(9)/h 

is an average over the measure of the broken symmetry group. 



4.1 A 2D example 

Before treating the three-dimensional problem of interest, it is instructive 
to consider a two-dimensional analogue, which is much simpler, but already 
contains the essential elements of the problem. To be specific, we consider a 
QPC whose width varies as 

D{z) = D* + z^/R, ze[-L/2,L/2] (12) 

along the wire (see Fig. |^). For a finite radius of curvature R, there is only 
a single unstable periodic orbit (plus harmonics), which moves up and down 
at the narrowest point of the neck. The monodromy matrix Mppo of the 
primitive periodic orbit is given by 

( dr' dr' \ 

dr dp 

dpi dpi I 

dr dp / 

with eigenvalues 



1 + D*/R D*{1 + D*/2R)/p\ 
2p/R 1 + D*/R J 



(13) 



e±x = 1 + D*/R ± ^{l + D*/Rf-l, 



(14) 
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Fig. 4. Point contact (a) and straight wire (b) as limiting cases of a nanowire. The 
point contact is characterized by the diameter D* and radius of curvature R of 
the neck. For the straight wire, R — > oo. A periodic orbit is sketched with a solid 
line, other orbits (dotted lines) are not periodic. The classical motion in the point 
contact (a) is in general chaotic, while the straight wire (b) possesses integrable 
classical motion 



2x being the Lyapunov exponent of the primitive periodic orbit. There is 
no continuous symmetry present (/ = 0), so the original Gutzwiller trace 
formula may be used to find 

^.o(i^)=^f:T^^?^. (15) 

^h^E^^^ |smh(nx)| 

In the limit i? — > oo, the Lyapunov exponent x ^ 0, and Eq. (^5|) diverges. 
In this limit, the wire acquires translational symmetry (/ = 1), and Eq. (^ 
may be used to find 

5gi{E) _ 2mD* ^ cos(2nfc£;D* - 7r/4) 

n— 1 ^ 

The classical motion is integrable in this limit. 

For large but finite radii of curvature, one can employ semiclassical per- 
turbation theory in R^^: 



(18) 



5g^,,,{E) = y ^= ^, (19) 



WD* J-L/2 

Ignoring the l/i?-dependence of the prefactor, one finds 

,^s„fn. _ C{^nk{E)LyRn) +tS{^nk{E)LyR7T) 

^nk{E)L^/Rn 
where C and S are Fresnel integrals. This leads to a DOS 
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Fig. 5. (a) Conductance G 
and (b) DOS giEp) for 2D 
nanocontacts with a ~ 5 
versus the contact diame- 
ter D* . g is normahzed to 
the area A of the region. 
Sohd curves: semiclassical 
results based on the inter- 
polation formula; crosses 
with error bars: numeri- 
cal results obtained by a 
recursive Green's function 
method [Q. Lower curves 
in (a) and (b): R ~ \f\ 
upper curves (offset verti- 
cally): i? 170Af 



where we have defined the function 

C[x, y) = cos{x)G(v) - sin{x)S{v). (20) 

The Gutzwiller formula ( p^ ) may be expanded in a Taylor series around 
i? = 0, while the perturbation formula ( pj|) gives a Laurent series around 
R — oo. Combining the two, an interpolation formula valid for arbitrary R 
can be constructed |lO|| : 



C ( 2nkED* ' 



nfi^kE ^1 |sinh(nx)| 



55int(i?)-^T27— r-TTTT-M (21) 



The crossover from integrable to chaotic behavior in Eq. ( pi| ) is controlled by 
the dimensionless parameter 

a{E) = L/^/X^, (22) 

where A^; — 2-K/kE is the de Broglie wavelength of an electron of energy 
E. We refer to a as the quantum chaos parameter: for a ^ 1 the DOS is 
indistinguishable from that of an integrable system, while for a 3> 1, the 
DOS is that of a chaotic system. 

In Fig. H, the DOS calculated from Eq. dH) plus the 2D Weyl term is 
compared to the result of a numerical solution of the Schrodinger equation. 
Remarkably, the semiclassical result is seen to be quantitatively accurate even 
in the extreme quantum limit D* ~ Af, -R ~ \f- 
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(2,1) 



(4.2) 



(3,1) 



(4,1) 




(5.2) 



(6.2) 



Fig. 6. Periodic orbits of an electron 
in the narrowest cross-section of the 
neck, labeled where v is the 

number of vertices and w the wind- 
ing number. The length of an orbit 
is = vD* sm<j>vw, where (pvw = 
nw/v is the angle of incidence at a 
vertex 




Fig. 7. DOS at the Fermi 
energy for axially sym- 
metric 3D nanowires of 
parabolic shape versus the 
Sharvin conductance Gs- 
The different curves repre- 
sent contacts with various 
ratios of D/Lq (indicated in 
the inset), where D is the 
asymptotic diameter of the 
constriction and Lq its ini- 
tial length 



4.2 3D nanowire with Eixial symmetry 

For an axially-symmetric three-dimensional nanocontact, the periodic orbits 
(see Fig. ^) occur in one-dimensional families which fit into the narrowest 
cross-section of the contact. This problem was first investigated by Balian and 
Bloch , who derived the axially-sy mm etric analogue of Eq. (p^ ) . We can 



follow the procedure outlined in Sec. 4.1 to derive an interpolation formula 



describing the crossover from a long nanowire (/ = 2) to a short nanocontact 
(/ = 1) 0: 



fvwLlwC {ksLyw - 3t!7r/2, a{E)^/v sin (py 



E 



where fyw = 1 + 0iv — 2w) counts the discrete symmetry of the orbit under 
time-reversal, the Lyapunov exponent Xvw is given by 



(24) 
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Fig. 8. Force oscillations SF versus the minimum contact diameter D*: dashed 
curve: \ima^o{SF}; solid curve: limc,-,oo{aSF}. The result for a 2> 1 is consistent 
with the WKB calculation shown in Fig. H, while the result for a (integrable 
hmit) agrees with the result M for a straight wire 



and the remaining terms are defined in the caption of Fig. |6|. Eq. ( |2 
plotted in Fig. |^. Note the rounding of the peaks in the DOS in short contacts 



IS 



5 Universal force oscillations 

The characteristic amplitude of the sawtooth structure in the cohesive force of 
a gold nanocontact was found to be of order InN, independent of the contact 
area or shape • To what extent is the amplitude of the force oscilla- 

tions universal? To calculate the force from Eq. we need to make some 
assumptions regarding how the shape of the contact scales under elongation. 
First, we assume that the deformation occurs primarily in the thinnest sec- 
tion, which implies D*^L « const. Second, we assume that R oc L^, which 
implies a = L/^JXpR ~ const. Inserting Eq. ( p3| ) into Eq. (||), and taking 
the derivative, we find [HOl: 



e_F jL^^ fy^ smjkpLy^ - 3vn/2 + 7r/4) 



w—1 v—2w 



SF ^ „ ^ y y A™ _ 3„7r/2). (26) 

Q<g;l w—1 v—2w 

^„ ep / 0.58621, a < 1, 
rms(5F = — X <^ ' _i ' (27) 

Xp [ 0.36208 a ^, a » 1. ^ ' 

From Figs. ^ and ^ one sees that a < 1 for a realistic geometry, implying 
that indeed rms 6F ~ ep/Xp. 
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Fig. 9. Artist's conception of 
a propagating Rayleigh insta- 
bility, from Ref. El 



A cylindrical body longer than its circumference is unstable to breakup 
under surface tension |^,^ (see Fig. |^). How then to explain the durabil- 
ity of long gold nanowires [c.f. Fig. ^b)], the thinnest of which have been 
shown ||2^ to be almost perfectly cylindrical in shape? Let us calculate the 
quantum corrections ||l^ to the classical stability analysis. Classically, only 
axially-symmetric deformations lead to instabilities. Any such deformation 
of a cylinder may be written as a Fourier series 



/oc 
dqb{q)e">', 
-OO 



(28) 



where R{z) is the radius of the cylinder at z and Rq the unperturbed radius. 
The constant volume constraint leads to the condition 



1 

b{0) = -— dq\b{qr. 
no Jo 

In the Weyl approximation, the energy cost of the deformation is 



(29) 



-fc^i?o + -fc^ ) &(0) 



+ / dq 



Hkjp / nkpRo Skp \ 2 
"15" [ 4 9~ ' 



For the unperturbed cylinder, Eq. (H) yields 



00 oc 



^ ^ TTh^ ^ V 

W—1 v—2w 



fvw^vw /IT o /n\ 
COS{kELy^a — 6V7T/2). 



(30) 



(31) 



The effect of the deformation may be treated with semiclassical perturbation 
theory: 



1 

LRt 



Jo 



dz R{z) e 



(32) 
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Fig. 10. Density of 
states giEp) of a 
cylindrical wire (upper 
diagram) and stability 
coefficient a (lower dia- 
gram) versus the radius 
Ro of the unperturbed 
wire. The wavevector 
of the perturbation is 
qRo ~ 1. Dashed curve: 
Weyl contribution to a 



where 

'^Syw (z) 



n 

Expanding Sg up to second order in b{q) gives 



dqb{q)e">'. (33) 



w—l v—2w 



6(0)(cos6'„tu - ksLy^o sin0y^) 



Ro 



dq\b{q)\'^ sine. 



(34) 



where 9vw{E) = ksLyw — 2,v'k/2. 

Combining Eqs. (^0|) and (|4|), and using the constraint (|2^), one finds 
that the change of the DOS is of second order in 6, and contributions with 
different q decouple. The energy integral (||) yields 

/•OO 

n[b]^Q[Q]+ dqa{q)\b{q)\^ +0{b^), (35) 
Jo 

where the stability coefficient a{q) depends implicitly on i?o and temperature. 
If a{q) is negative for any value of g, then Q decreases under the deformation 
and the wire is unstable. 

Fig. ^ shows the stability coefficient and DOS at the classical stability 
threshold gi?o = 1 as a function of Rq. The quantum correction destabilizes 
the wire where the DOS is sharply peaked; but what is more surprising, it 
stabilizes the wire in the intervening intervals. With these results, we can con- 
struct a stability diagram for the wire. For a given temperature, the stability 
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Fig. 11. Stability 
diagram for cylin- 
drical nanowires at 
two different tem- 
peratures. White 
areas are stable, grey 
unstable to small 
perturbations. The 
quantized electrical 
conductance values G 
of the stable configu- 
rations are indicated 
by bold numerals in 
(a), with Go = 2e^/h. 
Right vertical axis: 
corrected Sharvin 
conductance Gs- Dot- 
ted curve: stability 
criterion in the Weyl 
approximation 



problem is now determined by two dimensionless parameters: qR^ and kpRo- 
In Fig. |ll|, regions of instability, where a[q) < 0, are shaded grey, while stable 
regions are shown in white. Note that many of the white regions of stability 
persist all the way down to g = 0, indicating that an infinitely long wire is 
a true metastable state if its radius lies in one of the windows of stability. 
The quantized conductance values of the stable cylindrical configurations are 
indicated by bold numerals in Fig. |ll|(a). Our stability analysis is consistent 
with recent experimental results for alkali metal nanowires [^. 
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